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Abstract 

We study the consistency of parton distribution functions in the presence of target mass cor- 
rections (TMCs) at low Q"^. We review the standard operator product expansion derivation of 
TMCs in both x- and moment-space, and present the results in closed form for all unpolarized 
structure functions and their moments. To avoid the unphysical region at x > 1 in the standard 
TMC analysis, we propose an expansion of the target mass corrected structure functions order by 



^ ■ order in M^/Q^, and assess the convergence properties of the resulting forms numerically. 
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I. INTRODUCTION 



The application of the operator product expansion (OPE) to the phenomenological study 
of Quantum Chromodynamics (QCD) has been very successful in the determination of the 
quark and gluon substructure of the nucleon. The OPE allows the formal separation of cross 
sections for high-energy processes such as deep-inelastic scattering (DIS) into perturbatively 
calculable partonic cross sections and nonperturbative contributions parametrized by parton 
distribution functions (PDFs). The factorization of the cross section becomes especially clean 
in the Bjorken limit, where the energy v and four-momentum squared transferred to a 
nucleon with mass M both become infinite, with the ratio x = Q"^ /2Mv fixed. 

In the Bjorken limit the DIS process becomes dominated by scattering at light-cone space- 
time distances z^z^ ~ 0, with the expansion made in terms of products of singular and 
non-singular terms around the light-cone. The singularities are isolated in the perturbative 
Wilson coefficients, while the non-singular terms are all the possible operators allowed by 
the underlying quantum field theory. The coefficient of the operators of lowest twist (where 
twist is defined as the dimension minus the spin of the operator) contains the most singular 
terms. Operators in the expansion with higher twist are less singular, and at large are 
suppressed by powers of l/Q"^. 

While this framework has met with considerable success in describing data at high ^ 
and large final state hadron masses = M"^ + Q^{1 —x)/x, many recent high-precision 
experiments [1] have been performed at lower energies, with down to ^ 1 — 2 GeV^, where 
the use of the asymptotic Bjorken limit formalism is more questionable. In addition to the 
strong coupling constant becoming large, at low the higher twist power corrections, 
which describe nonperturbative multi-parton correlations, become increasingly important. 
Furthermore, even at leading twist, there are corrections arising from purely kinematic effects 
associated with finite values of Q'^ jv^ = 4M^,x^/(5^, usually termed target mass corrections 
(TMCs) [2-6]. 

To perform reliable perturbative QCD based analyses which include data in the low 
region, a careful treatment of the subleading corrections is essential, and global PDF 
analyses [7-11] have only recently begun to take such effects systematically into account. 
Studies of quark-hadron duality [12, 13] have also strongly suggested that data at low W 
can be described (to within ~ 10-15%) by leading twist parton distributions. A more basic 
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question, however, is whether one can consistently define leading twist parton distributions 
in the presence of TMCs, that can be valid at low over the entire range of x. 

The first analysis to tackle this question was by Georgi and Politzer (GP) [2], who pro- 
posed taking TMCs into account by defining distributions at low in terms of the Nacht- 
mann scaling variable, ^ [14, 15], 



This leads to a specific prescription for removing TMCs from measured structure functions 
that has been used extensively in the literature [5]. 

Unfortunately, problems with the standard TMC prescription were soon realized [16-21] 
in the behavior of the target mass corrected structure functions in the vicinity of x ~ 1. In 
particular, functions expressed in terms of ^ over the interval < ^ < 1 necessarily extend 
into the unphysical region between the elastic limit ^ = = ^{x = 1) and ,^ = 1 for any 
finite value of [2]. This not only violates the conservation of energy and momentum, but 
also makes structure functions nonzero at x = 1, at odds with the expectation that leading 
twist functions should vanish at the elastic point [12, 22]. 

As a possible remedy, De Rujula et al. [3, 4] noted that in the threshold region analyses 
of data should not be performed in terms of leading twist structure functions alone, without 
also incorporating the effects of higher twist operators. They argued that a nonuniformity 
in the limits as n — oo and — ?■ oo renders the entire approach untenable at very low W , 
when higher twists exceed ~ nM"^ /Q"^ for the n-th structure function moment. 

Attempts were also made by Tung and collaborators [19, 20] to phenomenologically re- 
move the threshold problem by utilizing an ansatz to smoothly merge the moments in the 
perturbative region at large with their correct threshold behavior in the — > oo limit, 
although such a prescription is not unique. Steffens and Melnitchouk [23] extended this 
approach by proposing threshold-dependent distributions which exactly satisfy threshold 
kinematics at all Q^, at the expense of sacrificing the universality of PDFs in the presence 



Other approaches based on coUinear factorization, starting with the seminal work of 
Ellis, Furmanski and Petronzio [2 1], avoid the inversion of moments by implementing TMCs 
directly in momentum space within the parton model [25-28]. These, too, however, suffer 
from prescription dependence [25-28], or do not extend to all orders in l/Q^ [24]. In addition. 




and /i = — . 



(1) 



of TMCs. 
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even though they do not invoke distributions at x > 1, all these formulations nevertheless 
retain the problem of nonvanishing structure functions at x = 1. 

Given the desire to maximally utilize the recent precision structure function measure- 
ments at large x [1, 7-11, 29], as well as those planned for the near future [30, 31], there is a 
pressing need to address the question of TMCs and the consistency of parton distributions 
with mass corrections at finite Q^. A more reliable treatment of the high-x region at mod- 
erate is important not only in providing a better understanding of the quark structure of 
the nucleon in the deep valence region [32, 33], it is also vital for constraining cross sections 
at collider energies through the evolution to lower x at higher values [34, 35]. 

In this paper we revisit the problem of kinematic thresholds and PDF definitions in the 
OPE approach to TMCs, elucidating its shortcomings, and proposing an alternative method 
that addresses some of the problems inherent in the standard TMC formulation. In Sec. II we 
review the standard TMC approach, outlining the OPE derivation of target mass corrected 
moments and their inversion to x-space. We demonstrate explicitly the conflict of the usual 
inversion procedure with energy-momentum conservation, and illustrate its consequences for 
the X dependence of the structure functions as well as their moments. We propose a new 
method to compute TMCs in Sec. Ill, based on inversion of the moments order by order 
in M^/Q^, without having to introduce the Nachtmann scaling variable ^, and study the 
convergence of the series numerically. The advantages and limitations of this method are 
summarized in Sec. IV. Further technical details of the TMC derivations of moments and 
structure functions are provided in Appendices A and B, respectively. 

II. TARGET MASS CORRECTIONS IN THE OPE 

In this section we begin by summarizing the basic formulas for inclusive cross sections 
and structure functions, before outlining the main steps in the derivation of TMCs from 
the operator product expansion. We present results for the complete set of leading twist 
moments of unpolarized structure functions, and discuss their inversion to obtain the x 
dependence at nonzero M^/Q^. 

In the one-boson exchange approximation, the differential cross section for a lepton scat- 
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tering from a nucleon target is given (in the target rest frame) by 

where VL is the scattered lepton sohd angle, and E and E' are the initial and final electron 
energies, respectively. The lepton tensors L^^ and the coefficients r] depend on the type of 
boson exchanged (7,7^, Z, or W^) [36]. Denoting the initial and final lepton momenta by 
k and k' , respectively, and the momentum transferred to the nucleon hj q = k — k', the 
hadronic tensor is given by the commutator of electroweak current operators J^, 

iyA'^ = i- J d^ze''i-'{N\[J>'{z),r{0)]\N) (3) 

p ■ q 2p ■ q p ■ q 

where Fj (i = 1 — 5) are the structure functions of the nucleon, usually expressed in terms 
of the variables x and = — g^, and we adopt the convention 60123 = 1 [•>()]. The structure 
functions Fi and F2 are accessible in charged lepton or neutrino scattering through a product 
of vector currents, while F3 requires the interference of vector and axial vector currents. The 
vector structure functions F4 and F5 are accessible in neutrino scattering but are suppressed 
by lepton masses, mf /M'^; for completeness, however, we include them in this analysis. The 
hadronic tensor Wfj_i, is also related to the imaginary part of the virtual forward Compton 
scattering amplitude, 

T"'' = 1 J d^ze"'-' {N\T{J^'{z)r{0))\N), (5) 

with 

W^^'^ = -discT'^^ (6) 

71 

In the following we derive expressions for the amplitude T'^^ and use Eq. (6) to extract the 
results for the structure functions. 
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A. Moments of structure functions 



The standard derivation of TMCs in the OPE in the twist-2 approximation begins with 
the Compton scattering amphtude T^^, which can in general be written as [2] 



k=l 

+ {g'.yq,. + g^.q'q,.) cf ) q,,--- g,,, ^ a,, w--^^'^, (7) 

where 

n™ = i2(-^y^^^i9---gp---PhApy (8) 

and {g - ■ - gp- ■ ■p}k,j represents the (symmetric) sum of {2k)\/[2^j\{2k — 2j)\] distinct prod- 
ucts of the form (7^n^»2 ■ ■ ■ gf^^2j-it^i2j ptJ-2j+i . . .ptJ-2k resulting from permutations of the indices 
/^i, ■ ■ ■ ) /^2fc- The Wilson coefficients Cf' are calculated perturbatively, while the factors 
are matrix elements of local twist-2 operators (9^i'"^2fe |^37j^ 

(iV| = A2kP^'---p^^^ - traces, (9) 

which parametrize the nonperturbative structure of the nucleon. In the case of the flavor 
singlet operator, for example, one has 

= iJ-f^'^'D^'^ ■■■Df'''^^ - traces, (10) 

where the braces {■ ■ ■ } denote symmetrization with respect to the indices /ii, ■ ■ ■ , fi2k- 
The Cornwall-Norton moments m}"^ of the structure functions Fi are defined by 

"1 

dxx''-^Fi{x,Q^) if z = 1,3, 4, 5 
Mt\Q')=l ' (11) 

[ dxx''-^Fi{x,Q^) ifz = 2,L, 
. Jo 

where the longitudinal structure function Fl is given by 

Fl = {1 + 4/ix2)F2 - 2xFi, (12) 
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with the corresponding coefficient function C£ = — C". A straightforward but tedious 
calculation gives for each of the moments [2, 38] 

MtW) ^ t (« ; {^Cr-' H- ^^,,J^,^_,f r-) A... (13a) 



j=0 
oo 



i=o 



M<»>(0^) = ± (" ^ ') (cr^- H- __4^_cr« ) A.«, (13c) 



i=o 



J ; V (n + 2j)(r2 + 2j-l) 



Mt' (Q^) = E /^^' ( " I ' j CT'' (13d) 



oo 



J ; ^(n + 2j)(n + 2j-l)'^2 + 4^^^ 



j=0 

-I^TTOJ3Ty'^r-)A.«, (13e) 
Mf '(Q^) ^ g (" ; ^) ^ (-^^^ CJ«^ . ic?«^) K,._, (131) 

where the binomial symbol (^) = a!/ [6! (a — 6)!]. Further details of the derivation of Eqs. (13) 
are given in Appendix A. Note that the expression for the M'f^ moment is the same as that 
in Ref. [38] once the differences between the corresponding operator definitions are taken 
into account [•>!)]. 

Up to this point the effects of the target mass on the structure function moments are rig- 
orously derived within the OPE formalism. To proceed beyond Eqs. (13) and determine the 
TMC effects on the x dependence of the structure functions themselves requires additional 
assumptions, which inevitably introduces some model dependence in the calculation, as we 
discuss next in the following section. 

B. Parton distributions with TMCs 

In the absence of color interactions, the matrix elements in Eq. (9) should not depend 
on any scale apart from the factorization scale. With this in mind, the products CP A^k in 
Eq. (7) can be written in terms of parton distribution functions /j as 

CfA^^ = [\yy''~'m, (14) 







where for ease of notation we suppress the dependence in CP and fi on the scale Q^, which 
arises from perturbative QCD corrections. The functions fi are defined such that in the 
massless hmit (/i — )■ 0) one has 

Ff, Ff , Ff , Ft\ Ff , Ff } = x/„ - A), /s, \f. ^h] , (15) 

where F-''°^ = hm^^o -^j is the massless limit of the physical structure function Fj. Note that 
our notation for the parton distribution functions fi differs from that in Refs. [2, 27], whose 
distributions effectively correspond to fi{x)/x. 

The functions fi can in principle be identified with the PDFs measured in deep-inelastic 
or other high-energy scattering processes. (For simplicity we omit the flavor dependence 
of the structure functions, including their electroweak couplings, which can be incorporated 
straightforwardly with the distributions fi.) Following the derivation of GP [2], the structure 
functions at finite can be inverted using the inverse Mellin transform, 

/ dnx-''M'f\Q'^) if 2 = 1,3,4,5 

2vr2 J_i^ 

F,ix,Q') = \ (16) 
/ rfnx-"+ii^^"^(g2) if i = 2,L. 

Using Eqs. (13) and (14), the x dependence of the structure functions can then be determined 
in terms of the functions fi, as outlined in Appendix B [2, 27], 

F.M^) = '^^^ ( ^^^) , (17b) 



dx^ Ve(l + Ai^ 



+ (1 + 4/ix )x ^.^ , , (17c) 



^3(x,Q^) = -x|-f-M)), (ird) 



dx \l + 

Fa[x,Q ) = — — pjrfiiO + I-IX — 



4(1 + /i^2) ^-v-^/ r + 

, ,V|1 . (He) 
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where the functions hi and gi are given by 

h.iO = f'dul^, (18) 
gM) = duhi{u). (19) 

Note that the expression for the -F4 structure function in Ref. [27] contains ".^/is" instead of 
"(75" in the second term of Eq. (17e). Equations (17) define the complete set of unpolarized 
structure functions in the standard treatment of TMCs in the OPE. As was noted aheady 
in Ref. [3], however, the standard resuhs lead to problems in the limit as x — ?■ 1, which we 
shall focus on in the remainder of this section. 



C. Consistency of the standard TMCs? 



When taking the moments of the calculated x-dependent structure functions in the 
presence of TMCs, one should for consistency recover the expressions for the moments in 
Eqs. (13). To be specific, we investigate this here for the F2 structure function, Eq. (17b), 
but the same arguments can be applied to all the other structure functions. From the 
definition of the moments in Eq. (11), the ra-th moment of F2 can be written as 



Mr(g 



(20a) 



d 



X 



dx V^(l + 

+ n{n - 1) / dxx'^~'^ 
Jo 



nx 



-1 3^^2(0 

e(i+/iaj.=o 



(20b) 



e(i+/ia 

where integration by parts has been performed twice. Changing variables from x to ^, and 
using the fact that the kinematic maximum value of ^ is given by .^o; the moment becomes 



+ nin — 



2vi92{io) 



n 



^2(^0 



1) E 



j=0 



n + j 



€0 



n+2j-2 



g2{0, 



(21) 



where we have also used dx/d^ = (1 + /i^^)/(l — fi^'^Y, together with the relation 



2\n+ 



n + j 



(22) 



= t diC^''-'92{0 + C diC^'=-^92{0- (23) 

^1 Jo J^o 



Now, consider the last term in Eq. (21) involving the integral of the function g2{0- From 
the definition of the parton distributions in Eq. (14), one can write 

1 r^o 
(n + 2j)(n + 2j 

However, because the function /2 (and hence its integrals as in Eqs. (18) and (19)) has no 
reason to vanish in the region < ^ < 1, the second term in Eq. (23) is in general nonzero. 
The same is true for the first three terms in Eq. (21), and as a consequence one does not 
recover exactly the original expression, Eq. (13b). 

On the other hand, if the parton distributions were to vanish in the region C,o < C, < ^, 
the moments would have to depend on ^q, 

C7A^{io)= f da-'m,i,) ^ ^MA= fd^^d^lMA ^ 0, (24) 
Jo "so Jo "so 

where we explicitly label the dependence of the functions fi on and .^q. This result suggests 

two immediate problems: (i) universal (process-independent) parton distributions would no 

longer exist at finite Q^; and (ii) the separation between short and long distances on the 

light-cone, as embodied in the OPE, would no longer be possible. 

If the condition that the structure functions vanish for ,^ > is not imposed, one is then 
faced with the prospect of energy-momentum not being conserved. In fact, if the upper limit 
of integration in Eq. (20a) were extended from x = ltox = l/(l — /i), the first three terms 
of Eq. (21) would be identically zero, and extending the integration in the fourth term to 
^ = 1, Eq. (13b) would be recovered. Consequently, the consistency of the GP prescription 
[2], and in most subsequent TMC treatments, requires the violation of energy-momentum 
conservation. It thus appears a general consequence of defining parton distributions at 
finite in the presence of TMCs that one must choose between two less than ideal options: 
either keeping a universal parton distribution and violating energy-momentum conservation, 
or conserving energy and momentum but working with process-dependent distributions. 

We can assess the numerical significance of the ^ > ^o region by evaluating the lowest 

(2) 

{n = 2) moment M2 using a simple form for the parton distribution, 

35 

xMx) = -y^(l-x)^ (25) 

chosen to approximately reproduce a typical valence quark distribution, normalized such 
that Jq dx f2{x) = 1. The moment shown in Fig. 1 is computed from Eq. (13b) for several 
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FIG. 1: n = 2 moments of the F2 structure function, illustrating the convergence of the series in 
Eq. (13b) for j = (dotted), j < 2 (dot-dash-dashed), j < 3 (dot-dashed) and j < 4 (dashed), 
compared with the standard TMC result from GP [2] using Eq. (17b) with the upper limit of 
integration ^max = ^0 (21) (dot-dot-dashed) and ^max = 1 (solid). 

values of j (from the leading term only, j = 0, up to the inclusion of the first four terms, 
j < 4), and is compared with directly integrating F2{x, Q^) over x from to 1 (or equivalently 
up to ^ = ^max = ^o); usiug Eqs. (17b). As noted above, this procedure does not recover 
the formal result for the moment, Eq. (13b), which is reflected in the nonconvergence of the 
moments with increasing j to the standard TMC result from GP [2] in Eq. (17b) (dot-dot- 
dashed curve in Fig. 1). On the other hand, if the missing integration range as expressed in 
the second term of Eq. (23) is kept, the convergence of the moments is recovered (solid curve 
in Fig. 1), although at the expense of effectively integrating beyond x = 1 (to ^^ax = 1)- 

The problem encountered here is at the core of the parton interpretation of the matrix 
elements An. The approach of GP attempts to maintain a partonic interpretation at finite 
by introducing a new scaling variable ^ [14, 15]. However, as shown in Eq. (21), this 
leads to inconsistencies in the extracted x dependence of the structure functions and their 
moments. 

A possible way to avoid the problematic ^ ~ ^0 region is, ironically, to not introduce 
the Nachtmann scaling variable ^ in the first place. This can be realized by performing the 
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inversion of the moments order by order in /x, rather than summing over all powers of // during 
the inversion. As we shall see in the next section, this allows us to work with universal twist- 
2 distribution functions, while simultaneously preserving energy-momentum conservation. 
The only drawback of this approach is that the region of x and where parton distributions 
can be formulated consistently in the presence of TMCs will be somewhat restricted. 



III. SERIES EXPANSION OF INVERTED MOMENTS 

In the course of inverting the moments to obtain the structure functions, the binomial 
theorem is used to perform the integration by absorbing combinatorial factors involving the 
integration variable n (see Appendix B). Instead of this standard procedure, in this section 
we describe how the moments can be inverted term by term by absorbing the combinatorial 
factor into derivatives, which gives rise to novel series expansions for each of the structure 
functions. We illustrate this procedure for the case of the F2 structure function, with the 
derivation of the other structure functions following similarly. 

For the j-th term in the series expansion for the F2 moment in Eq. (13a), , the 
contribution to the structure function is given by the inverse Mellin transform 

F2,A^.Q') = 77- / dnx'-^M^^{Q'). (26) 
27r« y_ioo 

Using integration by parts to write 



= (n + 2j)(n + 2j-l) f dyy''^^'-^ g^iv.Q 

Jo 



(27) 



the contribution to F2J can then be expressed in the form 



Next, we can observe that 



x-"-^^ =(-x)^+^-— (29) 



(n-2)! ' ' (9x2+J 

for all n on the imaginary axis except the origin, so that 

\2+j Q2+j rl rico 



-/-i^a^ I iv ^ y^'-^ ^ * (f) . (30) 
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Finally, making use of the 5-function representation in Eq. (B3) we arrive at the desired 
result, 



-X 



|2+i . Q2+j 



X 



This result can also be obtained by noting that, instead of Eq. (29), we can write 

l"- ^ J )■ n+2j-2 _ 2j ^ n+j 



y 



(n-2)! 

Substituting this into Eq. (28) then leads to 

1 1 



y 



y 



(31) 



(32) 



j!" io "^^^''^2(1/)^^ b'5(z/-a;)] 



(33) 



using again the 5- function representation (B3). Now, because the 5-function is a distribution, 
one defines its derivative (in analogy with integration by parts of regular functions) as 

dz S{z) (p'iz) = - I dz 5'{z) 0(^) (34) 



(35) 



for a given function 0. Applying this definition (2 + j) times to Eq. (33), we find 



F2,J{x,Q^) = ^^^ 



.1)2+3 



x^ I dy 



[y'^^92{y)] y^ S{y-x), 



j\ Jo dy^+'J 
which gives a result identical to that in Eq. (31). 

Either of these two methods may be applied to the moments of the other structure 
functions to obtain the complete expressions for the unpolarized TMC structure functions 
in terms of power series in M^/Q^, summing over all values of j. 



Fi(x,Q2) 
F2{x,Q^) 



j=0 

00 

x^J2fi- 



-xY 
j\ dx^ 



X ^ i;xfi{x) +jg2{x) 



^ [x'^92{x)] 



x' /i^- 



'-xV 



j=0 



j\ dx^ 



[x^^ 2 {xf2{x) - xfi{x) + 4:jg2{x))] 



' [x'^hsix)] 



F^{x,Q^) = xY,f^' 



j=0 



i=o 



j\ dx^ 



j\ dx^+'j 



1 



x^^ ^ [ jU - l)92{x) + -xMx) - jg^ix) 



x^^ ^ ( -392{x) + ^xh^{x) 



(36a) 
(36b) 
(36c) 
(36d) 
(36e) 
(36f) 
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Note that the Nachtmann variable ^ does not enter in Eqs. (36), and all the functions fi, Qi 
and hi are expressed as functions of x only. As required, the j = term in the expansion of 
Fi is simply the massless limit structure function, Ff^\ 

The advantage of this formulation is that it explicitly avoids the problems encountered 
with the consistency of the inversion in the GP approach discussed in Sec. II C. Indeed, 
direct integration of the structure functions in (36) leads to the correct expressions for the 
moments in Eqs. (13). 

To examine the convergence of the series in Eqs. (36), we show in Fig. 2(a) the first few 
terms in the expansion of F2, starting with the leading order, j = 0, term and up to the first 
five terms in the series, j < 5. For illustration, we use the simple massless limit function 
in Eq. (25), and compare the result with the standard TMC calculation from GP [2] in 
Eq. (17b). The results show that the convergence at = 1 GeV^ is fairly rapid for x < 0.5, 
with just the first two or three terms already giving a target mass corrected function that 
does not change noticeably with inclusion of higher orders. 

It is noteworthy that in this range one is already well within the nucleon resonance region, 
traditionally taken to be < 2 GeV, from which data are typically excluded in global PDF 
analyses. This can be more clearly seen in Fig. 2(b), where the structure function is shown 
as a function of W . The convergence of the TMCs is well under control down to values as 
low as ly ~ 1.3 GeV, just above the peak of the first resonance region dominated by the 
A(1232) resonance. At smaller W , or higher x, the higher order terms display oscillatory 
behavior as one approaches the nucleon elastic point, x = 1 (or W = M). For the particular 
form of /2 chosen in Eq. (25), x/2 ~ (1 — x)^, the first three terms in the series (j < 3) vanish 
as X —7- 1, while the contributions for j > 3 diverge at x = 1. The target mass corrected 
function from Eq. (17b) (labeled "GP" in Fig. 2) is finite at x = 1 and indeed extends into 
the unphysical region W < M. 

The large-x oscillatory behavior is significantly dampened by the time one reaches = 
5 GeV^, with the first three terms (j < 3) converging well up to x ~ 0.8, as shown in 
Fig. 3(a) for the ratio of F2 to the GP target mass corrected function, Eq. (17b). At this 

this corresponds to values of > 1.4 GeV, illustrated in Fig. 3(b), which again is well 
outside of the range where DIS data are typically used in global PDF analyses. Note that 
the vanishing of ratio of the leading order term, j = 0, to the full GP result as x — )■ 1 reflects 
the nonzero value of the GP TMC function at x > 1. 
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FIG. 2: Target mass corrected -F2 structure function at = 1 GeV^ from Eq. (36b), showing the 
convergence with increasing j, and compared with the standard TMC result from GP [2] using 
Eq. (17b), shown as a function of (a) Bjorken-x, and (b) the hadronic final state mass W . The 
arrows in (a) indicate the locations of the resonance region {W = 2 GeV) and the A resonance, 
while the arrow in (b) denotes the elastic limit, W = M. 
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FIG. 3: Ratio of target mass corrected F2 structure function at = 5 GeV^ from Eq. (36b) for 
various j {j = up to j < 5) to the standard TMC result from CP [2] using Eq. (17b), shown as 
a function of (a) Bjorken-x, and (b) the hadronic final state mass W. The arrows in (a) indicate 
the locations of the resonance region (W = 2 GeV) and the A resonance. 
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FIG. 4: Convergence of the series expansion for the target mass corrected F2 structure function 
for large values of j (j < 101, 102, 201, 202, 301 and 302), at = 2 GeV^. The arrows indicate the 
values of x at which the final state mass corresponds to the A resonance and to W = 1.25 GeV. 
Note the limited x range on the ordinate. 

One may ask whether the series converges over the entire range of x at finite if a 
sufficiently large number of terms is included in the sum over j. For the trial distribution 
(25) used here, Fig. 4 shows the result of summing up to ~ 100, 200 or 300 terms for 

= 2 GeV^. For x < 0.73 (or W > 1.27 GeV) fairly good convergence is observed, while 
for X > 0.74 (or < 1.25 GeV) the addition of a large number of terms is needed to push 
the convergence of the target mass corrected function to significantly higher x values. It is 
interesting to observe that the inclusion of additional odd or even terms in j gives alternating 
negative and positive divergent behaviors, respectively. The systematics of this convergence 
with and W will be discussed in more detail in Ref. [40]. What is clear, however, is that 
a severe limitation on the computation of TMC exists at very high x for small values of Q^, 
although this occurs deep in the resonance region at low W . 
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IV. CONCLUSION 



The problem of target mass corrections to deep-inelastic structure functions is almost as 
old as the theory of QCD itself. With the focus of most structure function analyses being on 
the perturbative region where subleading effects can be neglected, the study of TMCs 
remained largely dormant for several decades. The advent of new, high-precision data in the 
resonance-scaling transition region at high x and low has brought the problem of TMCs 
back to the fore, giving rise to greater urgency to the need for resolution of the remaining 
open issues with respect to their implementation. 

In this paper we have sought to illustrate the inherent problem with the standard TMC 
formulation, already evident in the pioneering work of Georgi and Politzer [2], in the treat- 
ment of the X ~ 1 region and the inversion of the structure functions from their moments. In 
particular, we have critically analyzed the definition of PDFs in the presence of TMCs, and 
discussed the consequences of the violation of energy and momentum conservation in the 
standard TMC analysis. Historically it has been argued [-3, 4] that the problem in the thresh- 
old region exists because at low the higher twist contributions cannot be neglected. We 
do not disagree that higher twists are essential for describing low energy structure function 
data; we believe, however, that one ought to maintain consistency of leading twist functions 
at any x, regardless of how large the higher twists may be at a given . 

We contend that the introduction of the Nachtmann variable ^, which appears naturally 
in the standard TMC implementation, does not lead to self-consistent parton distributions 
that are valid at all x. In fact, our analysis suggests that strictly speaking PDFs cannot be 
defined consistently at any finite when the mass of the target is incorporated. What is 
feasible, however, is to compute the x dependence of the TMC structure functions in terms 
of standard PDFs as a series whose convergence can be studied as a function of x and Q^. 

To this end, we have derived formulas for the entire set of unpolarized structure functions, 
as a series in M^/Q^, involving PDFs and their derivatives. The virtue of this approach is 
that the resulting TMC functions can be consistently inverted from their moments, with- 
out ever encountering unphysical regions of kinematics or violating energy and momentum 
conservation. Moreover, it allows us to systematically study the regions of x and where 
TMCs can be reliably applied. Using a simple trial function, we have illustrated the conver- 
gence of our scheme numerically for the case of the F2 structure function. Rapid convergence 
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is observed for most of the range of x, with the first two or three terms saturating the sum 
well into the nucleon resonance region. For values as low as 1 GeV^, we find that the 
convergence of the TMC series expansion is under control down to ~ 1.3 GeV, which is 
almost in the vicinity of the A resonance peak. 

At smaller W, or higher x for fixed Q^, rapid oscillations ensue as one approaches the 
elastic scattering limit, and beyond W ^ 1.3 GeV it becomes prohibitively difficult to 
tame these with a finite number of higher order terms. At = 2 GeV^, for example, 
even summing over ~ 300 terms allows for a smooth TMC function up to x ~ 0.73 (or W ~ 
1.27 GeV). Fortunately, such low W values are well outside the range typically encountered in 
perturbative QCD analyses of DIS data, and in practice will not pose any serious restrictions. 

Our results therefore lend greater support to global PDF fits which incorporate \ow-W 
data, currently down to = 3 GeV^ in some analyses [7-9], but with more ambitions 
plans to extend the range further into the traditional resonance region. Provided higher 
twist and other subleading corrections are tractable, our method for accounting for TMCs 
will introduce only minimum theoretical uncertainty into global analyses. A more detailed 
discussion of the utility of the present approach, as well as its application to spin-dependent 
structure functions, will be presented in a forthcoming publication [40]. 
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Appendix A: Derivation of structure function moments 



In this appendix we illustrate the derivation of moments of structure functions in the 
presence of TMCs, using as an example the Fi structure function (the Fi case contains 
some more general features that are not present in the F2 derivation discussed in Sees. II 
and III). The results for the other structure functions follow in a similar manner. 

We begin by finding Ti, the coefficient of —(7'"^ in Eq. (7), which has contributions from 
both the Cf' and Cf*^ terms. For the Cf' term, for fixed /c G N, j G {0, ■ ■ ■ , fc}, and each 
term of {g ■ ■ ■ gp- ■ 'P}k,j, a total of j q^^ factors will have their indices raised by j g^^^'- 
metric tensors. The result will then contract with j q^i^ to give a factor of {q^)^ ■ The 
remaining {2k — 2j) q^^ factors will contract with the {2k — 2j) p^' factors to give {p ■ qY^~'^K 
Since there are (2/c)!/[2''j!(2A; — 2j)!] terms m. {g ■ ■ ■ gp ■ ■ ■p}k,j, we find that 



k=l 



- g E(-i)- 2.,i'- m '"^ ■ iSr 

The g^^g''fj_^Q'^Cl^ term of Eq. (7) contributes to the coefficient of —g^^ due to the identity 
g^-^g'i^^g'^^'^'^ = g'^^ ■ For fixed A; G N and j G {0, ■ ■ ■ , /c}, we seek the terms oi {g ■ ■ ■ gp ■ ■ ■p}k,j 
which include a factor of gi'^i^^ ^jj-^g number of such terms is determined by the number of 
ways to distribute the indices /is, ■ ■ ■ , ^2k among {j — 1) g's and {2k — 2j) p's without creating 
duphcate products. From the (2/c)!/[2-'j(2fc — 2j)!] ways to distribute the indices /ii, ■ ■ ■ , fi2k 
over j g^s and {2k — 2j)p's without creating duplicates, relabeling indices j — t- j — 1 and 
k ^ k — 1 gives {2k — 2)!/[2''^-'^(j — l)\{2k — 2j)l] terms of {g - ■ ■ gp - ■ ■p}k,j which contain 
(^MiM2 fQj^ ]^ g j^^j ^ {1, ■ ■ ■ (Note that there are no terms of {g ■ ■ ■ gp - ■ ■p}kfl that 

contain g'^^'^'^.) As for the Cf' terms, we then find 

°° 22k 

E fe^.Q'C'f ) • ■ ■ q,,, A,, n---'^- 



k=l 

00 k 



k — 1 j — 1 

+ terms not involving g^^. (A2) 

Now, the term -te'^"''" g^^.q^q.^Cl" + {q^q" / Q'')q,,q,,Cf + {g^q^q,, + g'^Xq,,) Cf in 
Eq. (7) will not contribute to Ti, as the indices and u are "locked up" in such a way 
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that they cannot result in a factor of g'^'^ through contractions. We conclude, therefore, that 
the coefficient Ti of —g'^" in the expansion (7) is 

k — 1 j — 1 

(A3) 

Substituting p"^ = M^, = —Q^, and p ■ q = Q'^/2x into Eq. (A3), changing indices in each 
term to I = k — j and j, and rearranging, gives the result 

oo oo /ni A 1 

1=0 j=o \ J y 

Finally, using the identity §^ duj uj"'~™~^ = 27ri5„.m together with Eq. (6), we find that 



'c 

oo oo 



2 27rz In 
1 
2 



i=0 j=0 \ J / 

which leads to Eq. (13a). The results for the other moments (13) are derived in a similar 
manner. 



Appendix B: Structure function inversion 

In this section we illustrate the standard moment inversion procedure by presenting a 
detailed derivation for the case of the Fi structure function. The derivations for the other 
structure functions can be deduced straightforwardly from this example. 

We begin by denoting the series in Eq. (13a) involving and C^+^^ by mS"^(Q2) and 

m'^\Q'^), respectively. Using Eq. (14), we find for the m^i^ term, 

-I rioo -| rioo 
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A — n 



n+ j 
J 



if^yy- (Bi) 



From the (generalized) binomial theorem, it follows then that 

1 rioQ -| rioo pi -t 

dnx-''m^C\Q'^) = — dn dy x-"" y""-^ fi{y) 



where we have used the (5-function representation 

6i\nu) = — rfne™(i'^") = / rfnw". (B3) 

2vr 2711 J_i^ 



Using the relation 

a = root of u 

with u{y) = \n{y/[x{l — ^y^)]) and u'{y) = (1 + — /xy^)], the only root of u on [0, 1] 

corresponds to ^ = 2x/(l + p), which leads to 



For the mg""* term, using integration by parts with Eqs. (19) and (27), we can write 
— j dnx-mt\Q') = —j dn dyx-y-'g,{y)Y,ji^ ^ j {f^y^ ■ (B6) 
Next, from the relations 



we obtain 

— r dnx--m^^\Q') = ^f, r dn [' dy {n + 1) x'^y^ (B8) 



Finally, since (n + l)x ^ = —x'^{d/dx) x " ^, we arrive at the result for the m^^'' moment 



2ti J-,^ 27,, dx 7_i„ i„ " (1 - 



9a: Jq x{1 — iiy'^Y 



1 



dn 



(B9) 
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Combining Eqs. (B5) and (B9) then gives the final result for the inverted Fi structure 
function, 



The results for the other structure functions F2, - ■ ■ , -F5 in Eqs. (17) follow analogous deriva- 
tions. 
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